Proton-air cross section measurement with the ARGO-YBJ cosmic ray experiment 
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The proton-air cross section in the energy range 1-100 TeV has been measured by the ARGO- 
YBJ cosmic ray experiment. The analysis is based on the flux attenuation for different atmospheric 
depths (i.e. zenith angles) and exploits the detector capabilities of selecting the shower development 
stage by means of hit multiplicity, density and lateral profile measurements at ground. The effects 
of shower fluctuations, the contribution of heavier primaries and the uncertainties of the hadronic 
interaction models, have been taken into account. The results have been used to estimate the total 
proton-proton cross section at center of mass energies between 70 and 500 GeV, where no accelerator 
data are currently available. 
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I. INTRODUCTION 

In the last decades, high-energy proton-antiproton 
colliders allowed extending the study of hadronic 
interactions up to center of mass energy -/s = 1.8 TeV pj, 
0, Q ■ Moreover LHC experiments at CERN are expected 
to analyze proton-proton collisions and to measure total 
cross sections at -y/s = 14 TeV [41 . Anyway, at energies 
exceeding ^/s ^ 70 GeV, where total cross sections start 
increasing, the knowledge of hadronic interactions is 
limited by experimental uncertainties and/or lack of data 

The study of cosmic ray (CR) induced showers 
provides a unique opportunity to explore hadron 
interactions in an energy range which not only covers the 
LHC region, but also extends well beyond it, the natural 
beam of primary particles spanning up to extreme values 
Furthermore, some phase space regions cannot be 
easily accessed at colliders, where, as an example, most 
of the produced particles might escape undetected in 
the beam pipe in the case of soft, i.e. low momentum 
transfer, hadronic interactions. On the contrary, soft 
collisions play a major role in the development of CR- 
induced showers, carrying energy deep down in the 
atmosphere. In currently used interaction models, these 
collisions are treated by means of the Gribov-Regge 
Theory (CRT) Q via Pomerons exchange Q, while a 
perturbative QCD approach can be appropriate at most 
collider experiments. Finally, it must be considered that, 
in the case of CR experiments, the collision targets are 
light nuclei (i.e. N, O, etc., hereafter "air" nuclei) and 
nuclear effects cannot be neglected. These aspects make 
somehow complementary the two approaches and cause 
their comparison to be extremely interesting from the 
point of view of very high energy (astro)particle physics. 

In CR experiments, the proton-air (hereafter p-air) 
cross section can be measured in several ways depending 
on the energy range and detection technique. Suitable 
models can then be used to estimate the p-p cross section 
0. Because of the rapidly falhng CR energy spectrum, 
the use of ground based detectors is mandatory. This 
obviously implies that the method systematics are larger 
with respect to accelerator experiments, the interaction 
initial state being not under control. Moreover, as 
pointed out in [l3|, these experiments are actually 
sensitive only to the so-called production cross section 
defined as 

prod _ tot _ el _ qel /-, \ 

p—air p—air p—air p—air V / 

where (J^plair ^'^d (Tp_air the total and the elastic cross 
sections respectively, while (J^^lair refers to quasi-elastic 
processes in which the nucleus gets excited without 
secondary particle production. There is some ambiguity 
in the literature, since the quantity cr^^^j^, as defined in 
EqlU has also been referred to as inelastic or absorption 
cross section. Here we will follow the terminology 
introduced in [l3|. For the sake of simplicity, in the 



following we will refer to cTp^air simply as the p-air cross 
section or Gp-air- 

The p-air cross section is measured by evaluating 
the absorption of the primary proton flux penetrating 
the atmosphere. In principle, the distribution of the 
atmospheric depth, Xq, of the first interaction point of 
CR protons could give a measurement of the interaction 
length Xp-air (and then the cross section ap^air) m a- 
straightforward way. Unfortunately, this information 
cannot be currently accessed and indirect methods have 
to be used. Up to now three experimental approaches 
have been adopted, depending on the primary energy 
range and detection techniques. Their comprehensive 
description is out of the scope of this paper. Here we will 
shortly list them and we will then focus, in the following 
sections, on the technique adopted in this work. 

The first method dates back to 1960s and actually 
estimates the flux of protons that reach the detection 
level without having yet interacted producing an 
hadronic shower (see [ill. IH. [l3l| and references therein). 
This is made by measuring the flux of unaccompanied 
hadrons crossing high altitude calorimeters equipped 
with a suitable veto system. The flux is then compared 
with the measured primary cosmic ray flux at the top of 
the atmosphere. One of the intrinsic systematic errors of 
this approach comes from the uncertainty on the primary 
proton flux that is used in the estimate of the p-air 
cross section. By using this technique, one of the first 
evidences of the rising of the total p-p cross section 
with energy has been established [l3|- The reduction 
of the proton flux prevents the adoption of this method 
at high energies, where the experiments must rely on the 
observed extensive air showers (EAS) initiated by CR 
primaries at the top of the atmosphere. 

EAS arrays can evaluate the primary proton absorp- 
tion in the atmosphere by measuring the shower rates at 
different zenith angles (i.e. atmospheric depths) once the 
shower development stage and en ergy are constrained by 
observations at ground (see [H, [iQ, [13, US] and refer- 
ences therein). The measured number of muons, A^^, is 
generally used to evaluate the shower energy (and to se- 
lect proton- initiated events), while the number of charged 
particles at the detection level, N^, gives an estimate of 
the shower development stage. In this case an important 
role can be played by the shower-to-shower fluctuations, 
that might decrease the experimental capability of prop- 
erly measuring the p-air cross section [19i.i2Qi]. Since these 
effects depend on cosmic ray energy and mass composi- 
tion, on detector capabilities and location, and on actual 
analysis procedures, the sensitivity of each experiment 
has to be checked with detailed simulations. This tech- 
nique, traditionally known as the Ng — Nf^ method, has 
a further source of systematic uncertainty which comes 
from the appreciable model dependence of the predicted 
muon numbers pll . [2^ . 

At even larger energies, the shower longitudinal proflle 
can be directly measured with the air fluorescence 
technique (and/or using the lateral distribution of 
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Cerenkov light at ground). The atmospheric depth of 
the maximum shower development stage, X^^ax: gives 
information on the CR beam absorption, i.e. on the p- 
air cross section (see fl?', '23', '23] and references therein) . 
In this case, one of the largest systematics comes from 
the uncertanties on CR composition. This effect is also 
present in the previously mentioned approaches, but 
it becomes more and more important with increasing 
primary energy. In particular, the effect of CR primaries 
heavier than protons is negligible in the first approach, 
while it is small in the second one and a correction can be 
applied on the basis of CR composition measurements. 
At extremely high energies the CR composition is not 
sufficiently well known in order for that correction to 
be applied without introducing a sizeable systematic 
uncertainty. 

Once the p-air cross section is measured, theoretical 
models can be used to give the total p-p cross section. 
This conversion is essentially based on the Glauber 
theor y 12511 a nd has been discussed in several works 

In this paper we will report on the measurement of the 
production cross section between I-lOOTeV CR protons 
and air nuclei with the ARGO-YBJ experiment (see 
Secini. The analysis is based on the flux attenuation 
for different atmospheric depths, i.e. zenith angles, and 
exploits the high detector accuracy in reconstructing the 
shower properties at ground. This allows a different 
approach with respect to the standard — method 
used in EAS experiments that operate at larger energies 
and with lower shower front sampling capabilities. The 
results have also been used to estimate the total p- 
p cross section at center of mass energies between 70 
and 500 GeV, where no accelerator data are currently 
available [1]. The adopted strategy is discussed in Sec lIIIl 
while details on data selection and analysis procedure 
are given in SeclIVI Results are then given in SecfVl 
together with the discussion of the method systematics. 
Conclusions are reported in SeclVIl 



II. THE ARGO-YBJ EXPERIMENT 

Cosmic ray physics in the 1 — 1000 TeV energy range 
and very high energy (VHE) 7-ray astronomy are the 
main scientific goals of the ARGO-YBJ experiment, a 
full coverage extensive air shower array resulting from 
a collaboration between Chinese and Italian institutions 
[31I . [3^ . The detector is operating in the YBJ 
International Cosmic Ray Observatory, located in the 
village of Yangbajing, about 90 km north-west from 
Lhasa, in the Tibet region (People's Republic of China) 
at an altitude of 4300 m above sea level, corresponding 
to a vertical atmospheric depth of about 610g/cm^. 
The apparatus, extensively described in [3l|, [s^, is a 
single layer detector logically divided into 154 units 
called clusters (7.64 x 5.72 m^) each made by 12 Resistive 
Plate Chambers (RPCs) operated in streamer mode (see 



FigH]). Each RPC (1.25 x 2.80 m^) is read out by using 
10 pads (62 X 56 cm^), which are further divided into 8 
strips (62 X 7 cni^) providing a larger particle counting 
dynamic range '335. The signals coming from all the 
strips of a given pad are sent to the same channel of a 
multihit TDC. The whole system, in smooth data taking 
from mid 2006, provides a single hit (pad) time resolution 
at the level of 1 ns, which allows a complete and detailed 
three-dimensional reconstruction of the shower front with 
unprecedented space-time resolution (see Fig 12]). In order 
to fully investigate transient phoenomena with an even 
smaller energy threshold, an independent scaler mode 
acquisition system has also been put in operation l34i . 
Finally, a system for the RPC analog charge readout (S^ 
from larger pads, each one covering half a chamber (the 
so called big pads), is now being implemented. This will 
allow extending the detector operating range from about 
100 TeV up to PeV energies. 



A suitable calibration procedure has been developed 
in order to remove systematic time offsets in the read- 
out chain (due to front-end boards, TDC boards, cable 
lengths, etc.) among the 18480 different channels [36j . 
Many checks have been made in order to ensure the 
correctness of the procedure. One of these is made 
through the analysis of the position, size and shape of 
the reconstructed Moon and Sun shadows in the cosmic 
ray flux [13, HI. 



Due to the huge amount of data, event reconstruction 
is performed by using GRID-based object oriented 
technologies 



39| 



In the first event reconstruction 
step, the detected shower front is fitted to a plane. 
Actually, detailed simulations show that particles within 
several tens of meters from the shower core do form a 
curved front, which can be well approximated with a 
conical shape (40| . This is confirmed by experimental 
data and is suitably exploited in the second step 
of the event reconstruction. A reliable identification 
of the shower core position up to the edge, and 
slightly beyond, the active carpet is made through a 
Maximum Likelihood-based algorithm, using a NKG- 
like 41] lateral distribution function [44]. Space-time 
coordinates (positions and times of fired pads) are then 
fitted to a cone whose axis crosses the core position at 
ground [S]. As a result a very good angular resolution, 
fang IS obtained 38] . Obviously the actual value depends 
on the CR primary energy and mass (7, p. He, ...), and 
on the specific event selection criteria adopted in a given 
analysis. In our case, we get aang{E ~ 5TeV)~ 0.3° 
and (Tang{E > 10 TeV) < 0.2°, for primary protons. 
This is particularly important, since it implies that no 
significant systematics are introduced by the detector 
angular resolution in the zenith angle distributions used 
for the cross section measurement (see SeclV|. 
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FIG. 1: The ARGO-YBJ detector setup. The cluster (12 RPCs) is the basic detector and DAQ unit equipped with a Local 
Station for its read-out. The fuU-coverage central carpet is made by 130 clusters, the guard-ring by 24 clusters. The pad is the 
timing pixel and is further divided into 8 strips (see text). 
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FIG. 2: Two different views of a detected shower. The hit map at ground is given on the left, the color code representing the 
strip multiplicity of each fired pad. A space-time view of the detected shower front is given on the right. 



III. THE ANALYSIS STRATEGY 

For a primary energy interval and for a given distance 
Xdm between the detector and the shower maximum (see 
Figl3]), the frequency of showers as a function of zenith 
angle 9 is directly related to the probability distribution 
of the depth of the shower maximum P{Xmax), where 
Xmax = hosecB — Xdm and ho is the observation vertical 
depth. The shape of P{Xmax=Xo+Xrise) is given 
by the folding of the exponential dependence of the 
depth of the first interaction point Xq (i.e. Q-^o/>^p-air 
with Xp-air{g / cm?) ~ 2.41 X IQ'^ / ap-air{mb)), with 



the probability distribution of Xrisa which includes 
the fluctuations of the shower development up to its 
maximum. Then the effects of the limited experimental 
resolution have to be also taken into account. For 
sufficiently large X^ax values, P{Xmax) tends to be 
a simple exponential with a characteristic length A = 
kXp-air^ where k depends on hadronic interactions and 
on the shower development in the atmosphere and 
its fluctuations [1, [H, Hi, [ii]. The actual value of 
k must be evaluated with a full Monte Carlo (MC) 
simulation and it depends on the experimental approach, 
the primary energy and the detector response. It might 
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FIG. 3: Simplified scheme of the shower longitudinal development in the atmosphere (left panel) and of the Xmax distribution 
(right panel), together with the definitions of the quantities used in the text 



also depend on the features of the hadronic interaction 
model adopted in the simulation, this being a further 
source of systematic uncertainty. Anyway several studies 
showed that this dependence is generally small [^[islliH. 
Finally, the contribution of cosmic rays heavier than 
protons has to be estimated and minimized in order to 
get an unbiased proton- air cross section measurement. 

As discussed in SeclH experiments using the air 
fluorescence technique have direct access to P{Xmax), 
while EAS detectors measuring the particles at ground 
might sample it through the flux dependence on zenith 
angle, once Xdm (or the shower age) has been fixed or 
constrained, within the limits of detector capabilities. 
In these cases the combination of CR energy, detector 
vertical depth, angular and Xdm ranges, actually define 
the part of the Xmax distribution that can be accessed. 
A good performance is obtained if the exponential tail 
of P{Xmax) can be sampled. As it is schematically 
shown in FigO this requires selecting showers with the 
maximum development not far from the detection level 
(i.e. minimizing Xdm) and, obviously, exploring a zenith 
angle region as wide as possible. 

If this is not the flattening of the distribution 

might be observed, also due to shower fluctuations, 
resulting in a larger value of the parameter k and a lower 
sensitivity to Up-air- The measurement can be spoiled 
also by a very poor energy resolution. In this case both 
shower-to-shower fluctuations and the X^ise dependence 
on the primary energy might contribute to a worsening of 
the analysis performance. In some particular cases even 
the exponential behaviour can be lost ^2^. 

The ARGO-YBJ detector features and location (full 
coverage, angular resolution, fine granularity, small 
atmospheric depth, etc.), which ensure the capability of 
reconstructing the detected showers in a very detailed 
way, have been used to define the energy ranges and 
to constrain the shower ages. In particular, different 
hit (i.e. strip) multiplicity windows have been used 
to select showers corresponding to different primary 
energies, while the information on particle density, lateral 



profile and shower front extension has been used to 
constrain Xdm in the proper range. These features 
allowed the observation of the exponential falling of 
shower frequencies, through the sec9 distribution 

The fit to the angular distribution R[9) gives the slope 
value a, connected to the observed characteristic length 
A through the relation a = ho/ A, being 

i?(6i) = ^(6l)i?oe^"("""^-i' (2) 

where Rq = R{6=0), and the factor A{9) takes into 
account the geometrical acceptance of each angular bin. 

The same analysis chain is then applied to the 
simulated data sample (see next section). For each strip 
multiplicity (i.e. primary energy) interval, the fit to the 
sec6 distribution with the function of Eql^lgives the value 
of A^*-^ . The value of k , referring to each multiplicity bin, 
can then be evaluated as A: = K'^'^ / X^'Sair^ where Api'j^j^ 
is known, in the corresponding energy region, from the 
adopted interaction model. 

The experimental interaction length is obtained by 
correcting the observed characteristic length A'cR-air 
by the factor k determined on the basis of the MC 
simulation: A^^_^j^ = ^'cR-airl^- Such value will give 
the measured p-air interaction length (A^^^i^^j^), once the 
effects of heavier nuclei present in the primary cosmic 
ray flux have been taken into account. This has been 
made by evaluating the change of the slope a (see EqI2]) 
produced by the addition of the corresponding helium 
fraction to the proton primary flux in the MC simulation, 
the contribution of heavier nuclei being negligible. 

The proton-air production cross section is then 
obtained from the previously mentioned relation: 
ap^a^r (mb)= 2.41 X lOVAp!.';,^ (g/cm2), while several 
theoretical models can be used to get the corresponding 
total proton-proton cross section SecEl). 
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IV. DATA SELECTION AND MONTE CARLO 
SIMULATION 

The analysis was applied to a data sample of about 
6.5 X 10® events collected by the central part of the 
detector, i.e. the 130 adjoining clusters fully covering 
a surface of 78 x 74 (see FigH]). An inclusive trigger 
requiring at least 20 fired pads in a 420 ns time window 
was used. In order to have both a small contamination 
of external events (i.e. those events with the true core 
position outside the carpet but reconstructed inside) and 
an angular resolution better than 0.5°, only events with 
at least 500 fired strips were considered. Moreover, 
the analysis was restricted to events with reconstructed 
zenith angle 9 < 40°. This was made in order to avoid 
effects due to the possible zenith angle dependence of the 
analysis cuts, that might occur at 6* > 45° (see below). 

A suitable simulation chain was also set up to check 
the effects of the different analysis cuts and to estimate 
possible systematics. About 10® proton-initiated and 
2 X 10^ He-initiated showers, with the corresponding 
power law energy spectra between 300 GeV and 3000 TeV 
and zenith angle up to 45°, were produced with the 
CORSIKA code '46]. In order to have a better 
evaluation of systematics, we produced independent 
samples by using three different hadronic interaction 
models: QGSJET-I [13], QGSJET-II.03 gsf, and 
SIBYLL-2.1 [49, 53|. As reference values for the energy 
spectral indexes and t he p /He relative normalization we 
used those given in [5l| and resulting from a global 
fit of existing experimental data. A full simulation of 
the detector response, based on the GEANT package 
[S^ . was performed, including also the effects of time 
resolution, trigger logic, electronics noise, etc. Simulated 
data have been produced in the same format used for 
real ones and they have been analyzed by using the same 
reconstruction code. The reliability of the simulation 
procedure was successfully checked in several ways. As 
an example, in FigU there is a comparison between real 
data and the MC simulation, for two of the quantities 
used in the analysis As can be seen, the agreement is 
good both before and after applying the adopted event 
selection cuts (see below). Moreover, as expected, the 
agreement gets worse if a pure proton composition is used 
for the CR primary flux in the MC sample. 

In the analysis of real data, after a first selection 
based on the quality of the reconstruction procedure, a 
further rejection of external showers was performed with 
several additional cuts discussed in the following. The 
reconstructed core position, P^c, was required to be in 
a fiducial area given by the inner 11 x 8 RFC clusters 
(corresponding to a total surface of about 64 x 64 m^). 
This cut reduced the initial data set to about 45% of the 
events with reconstructed 9 < 40° and Nstrip > 500. The 
quantity i?7o was then introduced as the radius of the 
smallest circle (lying on the detector plane and centered 
in Pre) containing 70% of the fired strips. It was then 
required that the distance of Pre from the detector center 



plus i?7o should be less than 50 m. The aim of this cut 
was to select showers largely contained inside the detector 
area (thus very well reconstructed). As a result, the 
data sample was reduced to ~ 20% of the initial one. 
One further cut required the minimum average fired strip 
density, within a distance i?7o from the reconstructed 
core, to be 0.2strips/m^ in the shower plane. This 
density cut allowed the efficient rejection of external 
events, as shown in Fig. [5l left panel. The same purpose 
motivated the last selection cut (a compactness cut), 
requiring the i?7o radius to be at most 30 m. Monte Carlo 
simulations showed that this compactness cut is also 
related to the shower development stage, allowing the 
constraint of the value of Xdm (Fig. M right panel) and 
the selection of showers with their maximum lying deep 
in the atmosphere (i.e. the sampling of the exponential 
tail of the Xmax distribution). As already pointed out in 
Sec lIIIl this is an important point for the extraction of 
cross section data from the sec9 distribution. The last 
two cuts finally selected 12% of the events reconstructed 
with 9 < 40° and Nstrip > 500 in the initial data 
sample. The fractions of events surviving each analysis 
cut were checked to be consistent with the corresponding 
quantities for MC data. Indeed, once the contribution 
of CR primaries heavier than protons were considered, 
all the values were in agreement and independent of the 
adopted hadronic interaction model. 

The selected data sample was split into five different 
bins of strip multiplicity ANgtrip (see Tab|I]), starting 
from the threshold value of 500 fired strips on the 
whole central detector (out of the total 124800), in the 
trigger time window of 420 ns. As shown in TabUl each 
strip multiplicity bin corresponds to a different primary 
energy interval. In particular the averages of the Log(i?) 
distributions are reported (after the application of all 
the selection criteria) together with their root mean 
square values. The energy values E assigned in this 
way were checked to be equivalently given by the median 
energies E^q of the corresponding event samples, and to 
be independent of the hadronic interaction model used 
in the simulation. 

Simulations also show that, after applying all the 
analysis cuts, the contamination of external events 
misreconstructed as internal is less than 1% for all the 
energy bins. 

In order to avoid any bias in the cross section 
measurement, a check has been made that the event 
selection was independent of the zenith angle. As 
an example, the distributions of primary energies and 
Xdm of simulated events surviving the analysis selection 
criteria are shown in Fig|Sl for different zenith angle 
intervals. The distributions are independent of the 
zenith angle up to about 40°, thus showing that the 
experimental sensitivit y is not compromised by shower- 
to-shower fiuctuations |l9l . [20| . 
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FIG. 4: Comparison between simulated and real data before (left) and after (right) the application of selection cuts. The 
distributions of strip multiplicity are reported in the upper panels. In the lower ones the strip density near the shower core are 
shown, for events with the reconstructed core position inside the inner 11x8 RPC clusters (see text). The simulated sample 
include also helium initiated showers. As can be seen the agreement is good. Real and simulated distributions have been 
normalized, for comparison, to the same total number of events. 
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FIG. 5: Left panel: Correlation between the strip density near the reconstructed core and distance of the true shower core from 
the detector center. As can be seen, the density cut (dashed line) is able to efficiently reject misreconstructed events. Right 
panel: Correlation between Rjo and Xdm- The cut Rro < 30 m allows the selection of events with small Xdm (i.e. deeper 
shower maximum). These plots refer to simulated proton- induced showers processed as real data (see text). 
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FIG. 6: Energy and Xdm distributions (left and right panel, respectively) for proton-initiated MC events, corresponding to 
one of the strip multiplicity ranges, after the whole event selection procedure has been applied. As can be seen the effect of 
selection cuts does not introduce any substantial zenith-dependent bias. 



V. ANALYSIS RESULTS AND DISCUSSION 

The result of the application of the whole selection 
procedure to real data is given in Fig[7l where 
the experimental secO distributions for the five strip 
multiplicity ranges are shown, after correction for the 
geometrical acceptance of each angular bin. They clearly 
show the expected exponential behaviour, this being 
a further check that the detector capabilities and the 
adopted analysis cuts brought to a proper selection of 
events for the cross section measurement. A slight 
deviation is present, for the lowest energy sample only, at 
small sec9 values (therefore not included in the fit). This 
is interpreted as due to the effect of shower fluctuations 
(see Sec lIIip . the shower maximum being more distant 
from the detector for these events. Indeed, for the lowest 
energy sample {Xmax) ^ 390g/cm2, while (Xmax) ^ 
450g/cm^ for the other ones. 

The angular distributions obtained from the MC 
simulations are very similar to those obtained with 
real data and the same considerations can be applied. 
From these plots, the values of A*^*-^ are extracted 
by fitting them to an exponential function, with the 
slope parameter given by a = hf/'~^ / K^'^'-^ , where h^^'-^ 
is the nominal atmospheric vertical depth used in the 
simulation (see below). Such values, once divided by the 
values of \p'-airi give the parameter k for the different 
energies. In particular, for each of the five selected data 
samples, the averages of the corresponding distributions 
of Apl'jjj^ are taken, while their root mean square values 
have been used for the evaluation of the systematic errors. 

As mentioned before, we produced three independent 
MC samples with different hadronic interaction models, 
namely QGSJET-I, QGSJET-II.03, and SIBYLL-2.1, 



that are widely used for the simulation of EAS and 
have been shown to provide a consistent description of 
their properties [53| . In the 1-100 TeV primary energy 
range, their predictions for dp- air are quite similar, 
while there are larger differences for the estimates of 
the rate of diffraction processes and for the inelasticity 
of proton air interactions [H, H^l- This last quantity 
essentially gives the fraction of energy available for 
particle production [sl, [sl, [s^l and is particularly 
important for the longitudinal development of the shower 
in the atmosphere. Actually the two models that 
show the largest differences among them are QGSJET-I 
and SIBYLL-2.1, while QGSJET-II.03 gweii predictions 
closer to those of SIBYLL-2.1 (see for instance [1^). In 
order to have a comprehensive and conservative estimate 
of the systematic errors, we decided to consider all of 
them in the current analysis. 

The factors k evaluated with QGSJET-I, QGSJET- 
II.03 and SIBYLL-2.1 are reported in TabUl where both 
the statistical and the systematic errors are also shown. 
The first ones come from the fit procedure, while the 
second ones from the root mean square values of the 
corresponding A^l'^j^ distributions. As can be seen, there 
is a general agreement, with larger differences between 
QGSJET-I and the other two models, as expected. 
These results are also consistent with the lower average 
inelasticity of QGSJET-II.03 and SIBYLL-2.1, with 
respect to QGSJET-I iw the considered energy range (see 
below). By following the just mentioned conservative 
approach, the results have then been combined to get 
the estimated k factors and the associated statistical and 
systematic errors (last column in Tab|T|. 

As already discussed in Sec lIIIl the value of k results 
from the contribution of many effects that are difficult 
to disentangle and depend on hadronic interactions, pri- 
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N,„p = (500 -1000) 



20.38 / 8 
Constant 1 1.35 ± 0.01 
Slope -3.586 ± 0.036 



r^„jp = (1500 -2000) x'/nd1 24.8/13 
Constant 10.05 + 0.00 
Slope -4.61 + 0.03 
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FIG. 7; Experimental secO distributions for the five strip multiplicity ranges, after correction for the geometrical acceptance in 
each angular bin. 



Log(S/eV) 



kgasjET-iiA 



ks 



IBYLL-2.1 



k 



500 H- 1000 
1500 ~ 2000 
3000 ~ 4000 
5000 ^ 8000 
> 8000 



12.6 ± 0.3 
13.0 ±0.2 
13.3 ± 0.2 
13.6 ±0.2 
13.9 ±0.3 



1.98 ±0.06 ±0.05 
1.59 ±0.03 ±0.04 
1.69 ±0.05 ±0.03 
1.74 ±0.05 ±0.03 
2.04 ±0.06 ±0.05 



1.84 ±0.14 ±0.05 
1.75 ±0.12 ±0.04 
1.63 ±0.13 ±0.03 
1.97 ±0.17 ±0.04 
2.23 ±0.19 ±0.05 



1.87 ±0.08 ±0.04 
1.76 ±0.06 ±0.04 
1.72 ±0.05 ±0.03 
1.91 ±0.05 ±0.03 
2.01 ±0.05 ±0.05 



1.93 ±0.05 ±0.06 
1.63 ±0.03 ±0.08 
1.70 ±0.03 ±0.04 
1.84 ±0.03 ±0.10 
2.03 ±0.04 ±0.10 



TABLE I: Strip multiplicity intervals, corresponding proton primary energies, and correction factors k with their statistical 
and systematic uncertainties (see text). 



mary energy, shower-to-shower fluctuations, experiment 
location and capabilities, detection and analysis tech- 
nique. As shown in [HjIsHl, for ground array experiments 
a first contribution to k is given by: 



fco 



1 



1 - ((1 - i^„0''-'> 



(3) 



where 7 is the spectral index of the differential CR 
proton flux, Ki„ is the p-air inelasticity, and the symbol 
(...) indicates averaging over events. At our energies 
the inelasticity is expected in the range Kin — 0.53 
0.60 (depending on the assumed interaction model (s^ . 
Isgj). therefore implying fco ~ 1.25 ± 1.40. Further 
contributions to k are expected from a flattening of 
the sec0 distribution of selected events, produced by 
shower fluctuations and detector effects (see Sec lIII[) . thus 
making k larger. This can be parametrized by a factor 



kdet, being fc = fco fcdet [121,12 

As can be seen in TablH in our case fc ~ 1.6 ± 1.8 
(then we obtain kdet — 1-15 1.45 ), apart from the 
values obtained at the boundaries of the covered energy 
region. The larger value of k for the lowest energy bin is 
due to the smaller (Xmax) value, that produces a larger 
effect of shower-to-shower fluctuations (see Sec lIIip . An 
explanation for the relatively high value of k in the 
highest energy bin is given by the onset of saturation 
of the strip digital information used in the analysis. This 
makes wider the energy interval actually contributing 
to the considered multiplicity bin, therefore implying a 
larger effect of fluctuations, mainly in terms of Xrise, 
with a consequent loss of sensitivity. These two effects 
practically define the energy region in which the current 
analysis can be performed. Important improvements will 
be achieved with the use of the analog RFC readout (see 
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strip 

500 ~ 1000 
1500 ~ 2000 
3000 ^ 4000 
5000 ^ 8000 
> 8000 



V 

1.00 ±0.04 ±0.01 
1.00 ±0.03 ±0.01 
0.99 ±0.04 ±0.01 
0.98 ± 0.04 ± 0.03 
0.95 ±0.04 ±0.04 



CTp-air (mb) 

272 ± 13 ± 9 
295 ± 10 ± 14 
318 ±15 ±8 
322 ± 15 ± 20 
318 ± 15 ±21 



ap-p (mb) 
43 ± 3 ± 5 
48 ± 3 ± 6 
54 ± 4 ± 6 
56 ± 4 ± 7 
54 ± 4 ± 8 



TABLE II: Strip multiplicity intervals, correction factors for the helium contribution rj, and proton-air cross sections thus 
obtained. The values of the resulting total proton-proton cross sections estimates are also reported. The errors given first come 
from statistics, the second ones from the systematics. 



SecHH. 

A check was also made to verify that the detector 
angular resolution (see Sec|TT| is sufficiently good in 
order not to introduce any bias in the cross section 
measurement. Indeed the slopes a obtained by using 
the true zenith angles 9 or the reconstructed ones vi^ere 
in complete agreement, within the quoted errors, for all 
the five strip multiplicity MC samples. 

In the simulations the vertical profile of Linsley's 
standard atmosphere was used [i^ with the detection 
level set at the altitude of 4300 m a.s.l., corresponding 
to /iq^*^ = 606.7 g/cm^. A small source of uncertainty 
in the p-air cross section measurement (see Eq|2]) might 
be given by a systematic difference between h^'~' and 
ho together with its variations with time due to the 
change of atmospheric conditions. From the pressure 
data continuously taken at Yangbajing, we evaluated 
/ho = 0.988 ± 0.007, resulting in an impact on 
the cross section analysis at the level of 1%, which has 
then been taken into account for the cross section data 
presented here. 

As outlined in Sec lIIII the measured ^'cR-air value 
together with the k factor determined from the 
simulation, directly gives the experimental interaction 
length A^^_^j^ and consequently the production cross 
section acR-air- At this stage, the cross section has 
still to be corrected for the contribution of CR primaries 
heavier than protons. This correction has been estimated 
by evaluating the effect of the introduction of helium 
primaries in the simulated data on the shape of the 
sec9 distribution. Corrections for other primaries (i.e. 
CNO group, Fe, etc.) are negligible. Both the different 
absolute values and the energy dependences of proton 
and helium fluxes were considered, by taking as a 
reference the fits to the experimental data given in [5l| . 
As expected, the simulations showed a slight steepening 
of the sec9 distribution at the highest strip multiplicity 
intervals, thus changing the values of the corresponding 
cross section estimates. By applying this procedure we 
got the correction factors 77 to be applied to acR-air 
in order to obtain ap-air (see Tablll|. It has also been 
checked that the effect of primary particles heavier than 
protons on the shape of the sec9 distribution, is limited 
to few per cent not only due to the CR beam composition 
itself, but also because the analysis cuts actually select 



a proton-enriched sample. Indeed, the initial fraction of 
He-induced showers is further reduced by about a factor 
1.9, through the analysis flow, with respect to proton- 
induced events. 

The effect of the uncertainty of the primary CR 
composition has been estimated by applying the 
same procedure starting from other primary flux 
measurements, namely those from the JACEE[6y and 
RUNJ0B[61| experiments. We found small differences 
(below 4%) which were used for the evaluation of the 
systematic uncertainty on the factor 77 (see Tab HI)) . 

The resulting proton- air production cross sections, 
o'p-airj are summarized in TabHIl where both statistical 
and systematic errors are reported. These results are 
consistent with a previous ARGO-YBJ measurement 
done without using the strip information [6^ . and also 
with a preliminary estimate obtained from the analysis 
of a data set taken with only 42 RFC clusters during the 
detector installation [gl, [GJ] . 

The measured p-air production cross section is 
reported, in Fig|51 as a function of the primary proton 
energy. The results found by other experiments and 
the expectations given by some hadronic interaction 
models are also shown [iy, [gI, |6^. Even if data are 
in general agreement, it can be noticed some discrepancy 
between the various experimental approaches that might 
be due to the different method systematics (see SecH]). 
Furthermore, new results, namely those from ARGO- 
YBJ, EAS-TOPtlSil and HiRes systematically give 
cross section values that are slightly lower with respect to 
the more recent and comprehensive hadronic interaction 
models actually used in these analyses. The indication 
for lower cross section (and/or inelasticity) values is 
also consistent with what found from the analysis of 
several other EAS observables [63, [H, [6^. Indeed, 
as previously pointed out, the inelasticity distribution 
in p-air interactions plays an important role in the 
longitudinal development of atmospheric showers. In 
EAS-based cross section measurements these effects are 
taken into account by the k factors. Anyway residual 
systematic differences between data and expectations 
at the level of few per cent might result from the 
modeling and the effect of very low inelasticity collisions 
(i.e. Kin ^ 0.05) that cannot actually be observed, 
being due to diffractive interactions that weakly (or 
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FIG. 8; Proton-air production cross section measured by ARGO-YBJ and by other CR experiments [HI . [l^ . [isl . [l^ . [itI , [isl . 
m, m, [6^] together with the values given by several hadronic interaction models [4^ . |65| . IGo ]. Also shown are the predictions 
of two different calculations based on Glauber theory Isol ] applied to accelerator data as parametrized in [t^] • In this plot 
ARGO-YBJ points have been already corrected for the effects of CR primaries heavier than protons. 



do not) influence the EAS development [H, [H, [6^ . 
while are included in the cross section values given by 
the hadronic models in FiglS] [H, This kind of 

effects might be different in the case of cross section 
evaluations based on measurements of the single hadrons 
flux at ground [ill, 113, The low energy threshold 

of ARGO-YBJ (with respect to other EAS experiments) 
allows a direct comparison with the values given by 
this technique, showing a good agreement. This is 
particularly important since the systematics of the two 
measurement techniques are completely different. The 
agreement also extends to the predictions of different 
calculations based on the Glauber theory [2^, apphed 
to the measurements made at particle accelerators. As 
an example, the results of two of them [H, HO] are 
actually shown in Fig[51 starting from the accelerator 
data analysis performed in [t^I, where an asymptotic 
ln'^{s) dependence of hadronic cross sections is shown to 
be favoured. 

As outlined in SecUl Glauber theory can also be used 
to estimate the total proton- proton cross section Up-p, at 
energies not covered by accelerator data, by starting from 



the p-air production cross section (Tp^air as measured 
by CR experiments. These calculations were discussed 
in several papers 0, [13, [H, HE M, M, [13 . As 
pointed out by many authors [6|, an essential step in 
the conversion is the knowledge of the dependence of 
the slope, B, of the forward scattering amplitude for 
elastic p-p collisions on the center of mass energy ^/s, 
or on ap-p itself. For ARGO-YBJ data, we applied 
the conversion given in [l^ . that uses the somewhat 
model-independent Chou-Yang prescription [tH between 
B and the total cross section ap-p, this being also 
consistent with available measurements on B{s). The 
results were compared with what given by other models 
[2E[2l[ll[2l[la|, the larger differences being with [H 
and [23], that foresee slightly higher or lower ap-p values, 
respectively. Anyway in all the cases the differences, in 
our energy range, are below 10% and they have been 
added as a further contribution to the systematic error 
on the resulting ap-p. The results are summarized in 
Tabim and in FigO Also shown in the flgure are the 
measurements performed at accelerators [^ and the cr n_n 
values published by other CR experiments [l^, [IE [Ij] 
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FIG. 9: Total p-p cross section, o" p-p , obtained by ARGO-YBJ starting from a^-air together with the same quantitypublished 
by other CR experiments [l^, [2^.124]. The values of the total p-p cross section measured at accelerator experiments are also 
reported. The lines result from a global analysis of accelerator data [sl. ItoI. It^ . assuming an asymptotic behaviour of hadronic 
cross sections proportional to Zn^(s) or ln{s), as reported in the plot. 



starting from ap^air- As can be seen, the ARGO-YBJ 
data lie in an energy region not yet reached by_p-p 
colhders (and still unexplored by p-p experiments [l|, 0, 
1^). Our result is in agreement with the general trend of 
accelerator data showing the rise of the cross section with 
energy. In particular, it favours the asymptotic ln^{s) 
increase of total hadronic cross sections as obtained in 
[rot from a global analysis of accelerator data. 

VI. CONCLUSIONS 

The study of CR-induced showers provides a unique 
opportunity to explore hadron interactions in energy 
ranges and phase space regions currently not accessible 
at particle colliders. In this paper, the cosmic ray 
flux attenuation for different atmospheric depths (i.e. 
zenith angles) has been exploited in order to measure 
the proton-air cross section in a laboratory energy range 
between approximately 1 TeV and 100 TeV, with the 
ARGO-YBJ experiment. The information provided by 
a highly segmented full coverage detector, allowed the 



introduction of selection cuts able to both define the 
energy scale and constrain the shower development stage. 
The analysis results have been also used to estimate 
the total proton-proton cross section at center of mass 
energies between 70 and 500 GeV. 

The analysis was applied to a data sample of about 
6.5 X 10® events. A full simulation of the primary 
interaction, the shower development in the atmosphere, 
the particle detection at ground and the complete 
analysis procedure has also been setup and used in 
order to make consistency checks and evaluate the 
method systematics. The effects of shower fluctuations 
and the contribution of heavier primaries have been 
also considered. Moreover the use, in the simulations, 
of three different hadronic interaction models allowed 
a comprehensive and conservative estimate of the 
associated systematic uncertainties. 

The results for the proton-air production cross section 
have been compared with previous measurements and 
model predictions. A good agreement is present with 
other measurements made in the same energy region 
with completely different techniques, and also with 
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theoretical predictions based on particle accelerator data. 
Some systematic differences between more recent results 
and the values given by commonly adopted hadronic 
interaction models are present at the level of few per 
cent, the measured cross sections suggesting lower values 
and/or deeper shower development in the atmosphere. 

The total proton-proton cross section ap-p has then 
been inferred from the measured proton-air production 
cross section Op-air by using the Glauber theory. The 
resulting data lie in an energy region not yet reached by 
p-p colliders and so far unexplored by p-p experiments. 
Our result is consistent with the general trend of 
experimental data, favouring an asymptotic /n^(s) rise 
of the cross section. 

Further improvements in the analysis are expected 
from the use of the detailed information on the shower 
front (curvature, rise time, time width, etc.) that ARGO- 
YBJ is able to record with very high precision, and by 



the implementation of the analog RPC readout that will 
allow extending these studies to collisions with center-of- 
mass energies up to the TeV region. 
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